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Abstract 

We propose a Wigner quasiprobability distribution function for 
Hamiltonian systems in spaces of constant curvature — in this paper 
on hyperboloids — , which returns the correct marginals and has the 
covariance of the Shapiro functions under SO(D, 1) transformations. 
To the free systems obeying the Laplace-Beltrami equation on the 
hyperboloid, we add a conic-oscillator potential in the hyperbolic co- 
ordinate. As an example, we analyze the 1-dimensional case on a 
hyperbola branch, where this conic-oscillator is the Poschl- Teller po- 
tential. We present the analytical solutions and plot the computed 
results. The standard theory of quantum oscillators is regained in the 
contraction limit to the space of zero curvature. 
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1 Introduction 



In Hamiltonian systems which a have flat configuration space, among the 
phase space quasiprobability distribution functions, the Wigner function [p]], 
is the only one covariant under Euclidean translations of phase space 0, [|. 
The present paper and others that will follow it, aim to the construction 
of Wigner functions on configuration spaces that are conic surfaces, hyper- 
boloids and spheres, which transform under the Lorentz and rotation groups 
respectively, and which reproduce the traditional Wigner function when the 
conic contracts to the plane. Hamiltonian systems on conic manifolds have a 
natural kinetic energy given by the Laplace-Beltrami operator, and moreover, 
on these conies also a natural oscillator 'potential' can be proposed. In one 
dimension, this oscillator turns out to be one of the Poschl- Teller potentials 

I- 

In this paper, subtitled I, we propose a Wigner function on the D- 
dimensional hyperboloid H+, which generalizes the ordinary Wigner func- 
tion on flat phase space. It displays the correct marginals, and returns the 
traditional form of the Wigner function under Inonu- Wigner contraction to 
the zero-curvature limit. The elements and background for this assertion 
are contained in Section including the Shapiro solutions $' D) (x) to the 
Laplace-Beltrami equation 0. In Section |3| we present our proposed defi- 
nition of Wigner function on the hyperboloid, and verify the properties of 
marginality and the contraction limit to flat phase space. Covariance remains 
an issue because the Wigner function that we propose here follows from the 
covariance of the basis of wavefunctions, between the argument x and the 
index p, as if they were canonically conjugate variables. In this context, we 
re-examine the interpretation of momentum coordinates. 

In Section f| we exemplify the D-dimensional theory with a one- dimension- 
al sui generis oscillator on one branch of a hyperbola. This example may 
appear to be trivial, because the hyperbola is in most respects equivalent 
to a straight line. Nevertheless, the resulting Poschl- Teller potential is of 
particular interest because the wavefunctions are also the Clebsch-Gordan 
(Wigner coupling) coefficients for the three-dimensional Lorentz algebra, 
so(2, 1) = sp(2, 9R) = su(l, 1) 0. We display the Wigner functions of some 
Poschl- Teller wavefunctions; these have not been examined before. Finally, 
in Section || we recapitulate the aim and offer the present outlook of our 
program. 
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2 Elements of phase space and hyperboloids 



In his fundamental article |IJ], Wigner proposed a distribution function to 
represent on phase space the wavefunctions of pure and of mixed states in 
quantum systems. In this Section we recall the definition and properties 
that we shall generalize from flat to conic spaces, using the Laplace-Beltrami 
operator and the Shapiro functions. 

2.1 Wigner function on flat phase space 

In .D-dimensional flat configuration space x G 9? D , the generalized Dirac 
basis of plane waves solves the Helmholtz equation 

-A 0(x) = p 2 0(x), p (x) = exp(ip ■ x), pG^. (1) 

When we write p — +(p • p) 1 ^ 2 , n = p/p, and call p = pn the momentum or 
wavenumber vector, the functions p (x) represent plane waves in the direction 
of the unit vector n G <Sd-i in the (D— l)-dimensional sphere manifold. In 
the quantum model with natural units h = 1, p has units of inverse length; 
in the wave optical model, p is the wavenumber of light. 

The basis of plane wave functions ([!]) plays many roles: it provides the 
Fourier transform kernel which bridges the configuration and momentum 
realizations, it constitutes a basis for representations of the Euclidean group, 
and it serves for the construction of the 9ft 2£) -Wigner function of wavefields 
/(x), g(x) through the equivalent expressions 

Wx»(f,g\x,p) = — ^ / /V(x-iz)*e-* p - z s(x+!z) (2) 
= T^d I rf D z/(x-|z)*e +J P ^^) e - l P-( x +^)^ ( x+Iz) 

= 7Tw /rf D x'//x'7(xV^(x") 

x p( x') 5 D (x-i ( x'+x")) p( x")*. (3) 

This has the well-known properties of being sesquilinear in the functions, real 
for f = g, with the marginal projections J^d dpW = /(x)*g(x), J^d dxW = 
f(P)*9(P) (the tilde indicates ordinary Fourier transformation, F : f — /), 
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and covariant under translations in coordinate and momentum spaces 

T a : /(x) = /(x-a) => WV(T a :/, T a :^|x, p) = W^(/,£|x - a, p), (4) 
T b : /(x) = e ib x /(x) => V^ D (T b :J, T b : 5 |x, p) = <?|x, p - b), (5) 

F:/(x) = /(x) => F:<7 |x, p) = <? |p, -x). (6) 

The last intertwining by the Fourier transform was known when Garcia- 
Calderon and Moshinsky noticed that the Wigner function is covariant also 
under the larger group of Sp(2D, 9ft) linear canonical transformations of phase 
space 0. This is exceptional in the sense that the Heisenberg-Weyl algebra 
[whose generators are the phase space translations — and the unit 

that generates a commuting phase factor] has the outer automorphism group 
Sp(2Z},9ft). This accident does not occur for Lorentz algebras, so we should 
not expect similar covariances of the Wigner function under groups larger 
than SO(D,l). 

2.2 Laplace-Beltrami operator on the hyperboloid 

The purpose of this paper is to generalize the expression of the Wigner func- 
tion (0)-([3P with functions on a D-dimensional space x e 9ft D of constant 
curvature. This manifold can be seen in an 'ambient' space of D+l dimen- 
sions as a hyperboloid, with vectors x = (x ,x) e 9ft D+1 . 

Consider the upper sheet of the two-sheeted hyperboloid Ti® C 3ft D+1 of 
hyperbolic radius R > 0, 

, ,2 2 2 r>2 2 2,2, ,2 M\ 

\x\ = x Q — x = R , x = x 1 + x 2 + ■ — \-x D . (7) 

In this ambient Minkowski space, the isometry group is the Poincare group 
ISO(-D, 1)+, in place of the Euclidean group ISO(£>) + of flat space. The Lie 
algebra so(D, 1) has then the standard realization 

M j)k = Xjd Xk - x k d x . , M 0ifc = x d Xk + x k d xo , j, k = 1, 2, . . . , D. (8) 

The second-order Casimir operator, C, which is an invariant under the group 
SO(-D, 1)+, is (— R 2 times) the Laplace-Beltrami operator on TC 1 ?, namely 
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In the unitary irreducible representation spaces of the D-dimensional Lorentz 
group belonging to the most degenerate continuous series indicated by p 
[§], the operator (g) has a real lower-bound spectrum, as does ([!]). The 
wavefunctions of the free system on the hyperboloid are the solutions to the 
equation 



A LB /(x) = - 



D - l 

2R 

-- [0,oo 



+ p' 



/(X) 



x = -i(D-r 



A(A + D - 1] 
R 2 
ipR. 



/(x), 



(10) 



Any wavefield of a given wavenumber p is a solution of this equation. 



2.3 Shapiro functions 

A privileged basis for the solutions of the Laplace-Beltrami equation (|i0|) 
were given by Gel'fand, Graev and Shapiro |J in the form of D-dimensional 
plane waves of momentum p = p n, with positive wavenumber p and in the 
direction of a unit vector on the sphere n G iS D_1 , 



fxo-n- xX-K 5 - 1 )-^ 

I R ) 1 (11) 

(coshx - n • f sinh x )-^-i)-i pR ^ 



where functions f(x) on the hyperboloid x G H.+ (x 2 = R 2 ) will be denoted, 
according to convenience, by 

, ,/ x , x = + VR 2 + x 2 = R cosh x > R, / in \ 

fix) = /(x ,x) = /(x), x = R ^. nhx e ^ D ^ xe ^ G ^-r (12) 

The Shapiro functions (|TT|) are a Dirac basis for functions on the hyper- 
boloid, which are orthogonal and complete over x- and p-spaces: 

R r ePx 



[ IJ^^w = N ^(p)S D (p-vf), (13) 



(2tt) d Jxgsrd x 

1 f d D p 



r d D p (x) *&r» {x ' } = 5 D (x,x'), (14) 
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with the measure and Dirac 5 under J h d d D x = R J^d d D x./xo, 

T(ipR) 



N (D) {p) 
8 D (x,x') 



r(±(D-i) + ipR 



(pR) 



D-l 



<5 D (x-x'). 



(15) 
(16) 



In particular, N (1) (p) = 1, iV (2) (p) = coth(pR), and iV (3) (p) = 1. 

The Inonu- Wigner contraction limit of the Lorentz to the Euclidean group 
SO(D, — > ISO(D) + is the limit R — > oo in our expressions for vectors 
with x ~ R, x 2 <C i? 2 , and p = pn as before, i.e., 



lim $' D) (x) 



lim 

ii— >oo 



x — x ■ n\ 2 



R 



\{D-l)-ipR 



« lim [ 1 



x • n 



1? 



-ipR 



exp{i x ■ pj 



(17) 



Correspondingly, lim^oo N {D \p) = 1 and 5 (x, x') — > <5 (x — x'). 



2.4 Momentum space for the hyperboloid 

The Shapiro functions {$^ D) (x)} peK ,D in ([0]) serve as the integral transform 
kernel between functions of x on the hyperboloid, /(x), and conjugate func- 
tions of p, that has the interpretation of momentum or wavenumber space, 
and is indicated /(p). Using ( fL^) for x, p e 9?-°, one writes 

/(X » = A)*P° ,(X,/(P »- (19) 

This Shapiro transform has been used as a relativistic analogue of the Fourier 
transform (the physical context here, though, is not that of space-time rel- 
ativity, as we shall clarify below), and is a vector form of one of the two 
branches of the bilateral Mellin transform ||. Here the Shapiro transform 
replaces the traditional Fourier transform in the definition of a momentum 
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space p G 9? D , canonically conjugate with respect to this basis, to a con- 
figuration space of constant curvature. The corresponding Parseval relation 
is 

The manifold of momentum p = pn G (p G 3? + and n G 5' D_1 ), can 
be placed also in a (D + l)-dimensional 'ambient' space, where it occupies 
the cone w = (p, p) G V + - The momentum thus defined by the Shapiro 
functions has certain features however, which do not correspond to those of a 
standard relativistic momentum vector. If /(x) is a monochromatic wavefield 
with a definite value of p, this wavenumber will not change under SO(D, 1) 
translations of the hyperboloid ('boosts'), because it is the invariant value of 
the Casimir operator (fLOf) . Only the direction of momentum, n, can shift 
over the sphere; it will do so following the well-known Bargmann deformation 
of the circle [|K|], where the colatitude angle 'boosts' as tan ^0 i— > e _ ^tan|0 
for rapidity ( G 3?. Quotation marks are used for 'boost' because here we 
mean a translation in the hyperboloid, and not the well-known relativistic 
acceleration. 



2.5 Oscillators on conies 

The Laplace-Beltrami equation (|9|)- ([lOD provides the free fields (whose ener- 
gy is purely kinetic) on the hyperboloid. One introduces a potential energy 
term, as in the Schrodinger equations of quantum mechanics, by adding a 



function of position V(x) 0, 11, 12 



A LB + R 2 V(K)) / ( x) = R 2 Ef(x). (21) 

In quantum mechanics ji = m/h 2 , where m is the particle mass. For appli- 
cation in paraxial wave optics, we recall the interpretation where the extra 
term characterizes the refractive index anomaly of the medium, 

n(x) = n - z/(x), n o =n(0), \l (22) 

' v n i/(x) <-> V(x). v ' 
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A straightforward and useful generalization of the SO(-D)-isotropic har- 
monic oscillator potential from flat to conic D-dimensional configuration 
space is 



V(x) = i/iu; 2 ^ 2 ^- = \iiu 2 R 2 tanh 2 x = l^R 2 (1 - sech 2 ^), 



Xr 



(23) 



where x £ * s the hyperbolic angle coordinate defined in Eqs. (|T^). This 
is the Poschl- Teller 'secant-hyperbolic-squared' trough. 



3 Wigner function on the hyperboloid 

With the Shapiro basis of wavefunctions of the free system, we construct now 
our proposed Wigner function following the double-integral form in Eq. (^) 
for two wavefunctions, f(x) and g(x), by means of integrals on two hyper- 
boloids, \x'\ = R and \x"\ = R. 



3.1 Definition 

With the measures in Eqs. fll3|) and the Shapiro functions in ([□]), we define 
the Wigner function on the hyperboloid by 

W n (f,g\x,p) = — — -fi / — - / —j- f(x) g(x ) 

where A D (x; x', x") takes the place of the Dirac delta 6 D (x— ~(x'+x")) on flat 
space, Eq. (|]), and which will be detailed below. 



The crucial property that we must require of this 'binding-/!' in ( |24] ) is 
that it should guarantee that x be the midpoint of the geodesic between x' 
and x", so that all three points lie on the hyperboloid We achieve this 
in the following way ||13|| : given x G 3l D+1 in the upper sheet of a two-sheeted 
hyperboloid, we build any y G 3? D+1 on a one-sheeted hyperboloid HP of the 
same radius R, such that it be Minkowski-orthogonal to x, 

y = (yo, y), \y\ 2 = yl - y 2 = -R 2 , x y - x • y = o. (25) 



8 



Then, we can express x' and x" as vectors obtained from x and y as follows: 

x' = x cosh \t — y sinh ~r, x" = x cosh ~r + ?/ sinh ~t, (26) 

where x', x" G 7i+ for all r G 9ft. Also, it is easy to show that 

XqXq — x' ■ x" = R 2 cosh r, xox' — x • x' = xqx'q — x ■ x" = R 2 cosh |r, (27) 

i.e., the geodesic distance between x' and x" is i?r, while x is at |i?r from 
both x' and x". The arguments x' and x" in the expression fl2~4] ) thus emulate 
the arguments x ± in (^|) with the parameter |r. 

Using ([R]) and the parameter r in Q2"7|), we propose the binding- A in ([2~4]) 
to be 



This will yield the correct marginals (to be seen below) due to its properties 

A D (x;x',x') = — <5 D (x-x'), R f — A D (x; x', x") = 1. (29) 
it ixe» D Xo 

3.2 Integral forms 

The 2.0-fold integral form of the Wigner function in ([24]) contains Dirac 5's; it 
can therefore be brought to a (D + l)-fold integral noting that the definition 
of y G HP leaves the freedom of rotating y around x on a sphere S D ~ X . 
When we change variables from x' and x" to x and y according to (p6|), we 
reduce the integration to y and r while keeping Minkowski-orthogonality. 
The proposed Wigner function (E¥) then becomes 



R 2 f°° f 
W H (f,g\x,p) = 77TT5 (sinh r^dr / „ d D y S(x y - x ■ y) 

(27l) n JO JyeH D 

x f(x cosh \t— y sinh ~t) #(x cosh hr+y sinh |r) (30) 
x$j, D) (x cosh |r— y sinh |r)$| ) I?) (x cosh |r+ysinh |r) . 

The Dirac 5 remaining in (^) can be used to find a third alternative 
form of the Wigner function. This is obtained with the parametrization of 
the ambient-space vectors given by 

x = (x , x) = i?(coshx,£sinh x), V = (Vo, y) = R(smh u,rf cosh u), (31) 
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where £ and 17 are unit vectors on the sphere S D 1 and 9?o~. The Dirac 

5 in Eq. (|30|) is then 



L cosh 

£(xo2/o — x ' y) — — n -, — S(u — Q), with tanh Q = £ • ?/ tanh x- (32) 

it" 1 cosh x 

The differential oPy of the integral in y e becomes .R' D (cosha;) D ~ 1 x 
da; d D ~ l r], so the Wigner function ( pip becomes a .D-fold integral with the 
structure of (|2T), to., 



d-1 , /• iyi 



D 



{2k) u Jo Js - 1 cosh% 

x f(x cosh \t— y sinh |r) cosh |r+y sinh |r) (33) 

x $ jf' (s cosh |r— y sinh ir)$[ > D) (x cosh |r+y sinh |r) , 



with 



,/ f/lsmliO ,/cosh <>) = ^^ tanh ^ g) , (34) 

1 — (£ • r] tanh x) 2 



3.3 Marginal projections 

The marginal projections obtained by integrating the proposed Wigner func- 



tion (24) over momentum and configuration space should yield, respectively, 
f(x)*g(x) and /(p)*fif(p) as defined in fll8f)-fll9D . The two marginals follow 
from the orthogonality and completeness relations of the Shapiro functions, 
Eqs. (HK0). 

The integration of the Wigner function over §l D momentum space with 
the measure l/N (D) (p) in (|15D, is 



= R2 1 f*r d°^_ f 3 j* ff) a d ( x , xll) 5 d {xI j 

r d D x.' 

= R — i -f(x')*g(x')A D (x;x',x') 
A'eK° x 

= f(x)*g(x), (35) 
where we used (0), (0), and the first property of the binding- A in (|2 
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Similarly, the integration over configuration space with the measure 
R/x in (0), is 

M H (f,g\p) = Rf jL?— W H (f,g\x,p) 

= f<P)*9(P), (36) 

where we used the second property of the binding- A in (|29|) and the Shapiro 
transform (|18|). 

Finally, integrating over the whole of phase space with the appropriate 
measures in the Parseval relation fl20|), we have the total probability 



R 



f /^M H (f,g\x) = (f,g) H o= [ -^-M n (f,g\p). (37) 



3.4 Covariance under rotations and conic translations 

Under rotations R G SO(D), wavefunctions f{Xo,x.) transform through 

T(R):/ ( x) = /(a;o,R- 1 x). (38) 

In particular, the basis of Shapiro functions (Til"| ) transforms as 

T(R) : $; D n »(a:o,x) = $£>(z , R^x) = $^ n (x ,x). (39) 

Applying rotations T(R) to the wavefields / and g in the Wigner function 
(24), we next change variables to x' = Rx' and x" = Rx" (the ambient 
Xo-components behave as scalars), then use (|3~9"D for R _1 , noting that the 
binding- A in ( pip is invariant, A D (x;x',x") = A D (x;x',x") for x = R *x, 
and so are the measures d D x. r = cPx'. It thus follows that the Wigner 
function ( pi] ) is covariant under rotations, fulfilling 

W H (T(R):f,T(R):g\x,p) = W H (f,g\R-\R- l p). (40) 

Now consider translations by ( ('boosts' of rapidity () B m (£) G SO(D, ly + 
in the direction of unit m G S " 1 , which transform the ambient space vectors 
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preserving the constant-curvature subspaces x E for each radius R > 0. 
We denote by xy m and xi m the projections of x parallel and perpendicular 
to the direction of m, so that x = X|| m + x im . Then, wavef unctions on the 
hyperboloid transform as 

/ Xo \ I Xo cosh ( — m • x sinh £ \ 
T(B m (C)) : / X|| m =/ X|| m cosh C-£ m sinh ( • ( 41 ) 

\ x ±m/ V X ±m J 

When this transformation is applied to the plane-wave basis of Shapiro func- 
tions, their directions n on the sphere change, and they acquire a multiplier 
factor: 

T(B m (C)) : ^(i 0l x) = (coshC+m-n sinh C )-|P-D-*R (xQ) x)> (42 ) 

where the components of n' e S 0-1 that are orthogonal and parallel to m 
are 

n z _ n ±m ^ _ n ■ m cosh ( + sinh ( 

~ Lm cosh ( + n ■ m sinh ( ' " m cosh ( + n • m sinh ( ' 

within the same SO(D, 1)+ irreducible representation characterized by the 
invariant wavenumber, p. If the angle from m to n is 0, it will transform 
through the well-known Bargmann SO(2,l) map of the circle. 
The expression in the multiplier factor of Eq. 



/i(m, C; n) = cosh^+m-n sinh^, (44) 

is, not coincidentally, p'/p — if the (D+l)-vector w = (p, p) G V + were 
allowed to transform as a 'lightlike' vector in relativity, i.e., without being 
constrained to its p-sphere. Under the inner product (0), the SO(D, 1) 
boost with the multiplier (|42|) is unitary nonetheless, because the measure in 
p-space is d D p = p D dpd N ~ l n, and while p is invariant, from fl4"5| ) it follows 
that d D ~ l xi = /i(m, £; n)^^ 1 <i D ~ 1 n / . This cancels the absolute square of 
the multiplier ( f4"4] ) in the Shapiro functions fl4"T|). This type of covariance 
modulo a multiplier function is determined by the Shapiro function basis; we 



may call Eqs. ( |4T1) and ( |42D the Shapiro covariance between the conjugate 
transformations in x and p. 

When the wavefields in the Wigner function (|24l) are translated within the 



hyperboloid by (f41"D, the ambient vectors x are multiplied by a (D+l) x (D+l) 
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('boost') matrix B m (C), to x' ^ x' = B m (C)~V andx" i-> x" = B m (C)~V. 
Under this transformation, the measures are again invariant, d D x.'/x' = 
d D x.'/x' Q , etc., and so is A D (x; x', x") = A D (x; x',x"); hence, x = B m (^) _1 x 
will appear in the first argument of the transformed Wigner function. But 
the corresponding transformation of p in each of the two Shapiro functions, 
Eqs. ([12|) and (|44]), yields a multiplier factor. The imaginary exponents of 
/u(m, (, n) cancel, and there remains a positive net multiplier factor: 



W n (T[B m (()]:f, T[B m (Q}:g 



x,pn 



= (//(m,C,n)) ^Vw^lBK)- 1 *, pB^Q-'-.n), (45) 

where n' = B m (C) _1 :n is given by ([43|). Note that in the D = 1-dimensional 
case, the multiplier factor is 1. 

Covariance of the Wigner function is usually understood in the simple 
form it has under rotations, as given by fllCf ). Under these transformations, 
the hyperboloid in the ambient x-space rotates on its axis, and in the mo- 
mentum plane the circles n of all radii p rotate in sinchrony. Translations 
within the hyperboloid fl4li| ) on the other hand, deform the ambient and 
projected space vectors, x and x, through (f4~T|); momentum space is concur- 
rently squeezed in the direction of the translation so that its points move on 
constant-p circles and with a common Bargmann deformation of the angle. 
Since areas are not conserved in momentum space, a multiplier function of 
p is necessary for the Wigner function to ensure the total conservation of 
probability 



3.5 Contraction limit 

We now show that, when f(x) and g(x) are significantly different from zero 
only within a small, essentially flat patch of the hyperboloid, the definition 
of the Wigner function in Eq. (^) reduces to the standard Wigner function 
for flat space in Eq. ([|). In (|33D , the integrand for rj [recall Eqs. (|3lD and 
(|34"D1 will be significant only when 9ft D norms of the vectors fulfill 

i u i _l • i i i ^ d f |x| cosher < R sinhx < 1 , Aax 

xcosh ±r ± y smh <C R =>• < . , i „ . , , 46 

1 2J2I | | y | I r < _R =>. sm h r < 1 v ; 

xnR(l,xO, V~ K(\Z-v-v)- (47) 
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Also, using the limit in fllTD , and approximating sinhr ps r and cosh |r ps 
coshx ps coshcu ps 1, the Wigner function in Eq. (|33| ) reduces to 

HW, 9 |*,P) = r^drj^d^ (48) 

x /(x ,x— \Rrrf)* exp(-iRTT) ■ p) g(x ,x+^RTTj). 
Finally, changing variables to z = i?r 77 and noticing that 

/ d D z--- = R D r^-Ur f d - 1 !,--- (49) 

we see that PS|) reduces to (|2|). 

3.6 Special case of one dimension 

In the case D — 1, the Wigner function (^) actually coincides in form with 
the corresponding one-dimensional standard flat space form (0), as we now 
proceed to show. 

First, notice that the unit vectors n and £ in Eq. (|TTD are now the unit 
scalars n, £ — ±1, and that the Shapiro functions become simple exponen- 
tials: 

(R£ cosh X , R£ sinh x ) = [exp(-n£x)]~ ipR = exp(in£ P R X ). (50) 

We can now let np i— > p and £x i— > x with p, x e (— oo, oo), and recognize 
that floPP is a 1:1 function of only one variable, x\ G 3?, 

&p\xi) = exp(ixpR), xi = i?sinhx. (51) 

The arguments x = (xo, x\) of the functions / and g in ( p3[ ) then simplify, in 
components, to 

x =( x o\ x cosh i r±y sinh i r = R ( cos Mx ± M \ 

V^i/ 2 2 V smh(x ± ir/r) / 

For short, we indicate /(i? cosh x, R sinh x) = f(x)- The unit vector 77 in 
(|33l)-(|34T) also becomes a unit scalar, 77 — ±1, and the integral extends over 
y = r]i?(sinhx, coshx), and = tjx [see Eq. (|32|)1. Finally, the integral over 
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rj reduces to a sum over 77 = ±1, and for r r]r E (—00, 00), the Wigner 
function ( |33D becomes 



W H (f,g\x,p) = — drfix-lrTe-^gix+lr) (53) 

1 roc „ 

= — / dvf(p-\vye +iRv *g{p + \v). (54) 

The last expression is the usual flat-space Wigner function in terms of the 
conjugate wavefunctions on momentum space. Finally, note that for D = 1, 
the net multiplier which appears under 'boost' transformations in (|45|) is 
unity, so standard and Shapiro covariances coincide. 

4 Example: oscillator on the hyperbola 

We consider the open one- dimensional space which is the upper branch of a 
hyperbola of fixed radius R > 0, 

U\ = {(x ,xi) G I xl -x\ = R 2 }, (55) 

parametrized as usual by the hyperbolic angle x ^ 

4.1 Laplace-Beltrami operator and the oscillator 

When the potential V(x) is a constant (corresponding to a homogeneous 
optical medium), the D = 1 Schrodinger equation (pl~l) is the free wave 
equation, and its Shapiro solutions are simply the oscillating exponentials 
(|5lD, with energy E = p 2 /2fi > 0. 

Since the Laplace-Beltrami operator on the hyperbola (|5"5]) is A L 



R 2 d 2 /dx 2 , the Schrodinger equation for the conic oscillator (|23|) is 

— 1 d 2 



2/i dx 



2 



-R 2 E sech 2 x)f(x) = R 2 (E-E )f( X ), E = \^ 2 R 2 , (56) 



V( X ) = \^ 2 R 2% -\ = ^s(s + l)(l-sech 2 x), (57) 

Xn 



\ + J + i (58) 
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The bound solutions are 



2~(s-n) 



[s-n)T(2s-n + l) 



r(s-n + l) V 
xsech s_n x 2^1 



n! 

— n, 2s — n + 1 
s — n + 1 



1 — tanhx" 



(59) 



-n) n\ 



\ 7rr(2s-n+i; 



r( s -n+I)(2sech X ) s -"Q-" +1 / 2 (tanh X ), 



where n is a nonnegative integer bounded by s + 1, and C°(£) are the Gegen- 
bauer (or ultraspherical) polynomials [15| for a > — |. The corresponding 
quantized values of the energy are quadratic in n, and counted from the 
lowest level up by 



K 



^ 2 R 2 



2/ii? 2 



n — s 



n 



0,1,2, 



< s + 1. 



(60) 



As a check on our concepts, we verify that the contraction limit R — > oo 
of this system, when the radius of the hyperbola grows without bound, is 
the harmonic oscillator on flat space. Since the coefficient s correspondingly 
grows as s ~ fiuR 2 , the linear-quadratic spectrum of energies in Eq. (|6C|) 
becomes the linear spectrum of the quantum harmonic oscillator E n = u(n + 
|). To implement this limit on the wavefunctions (|59"D , it is convenient to 
use the following forms for the Gegenbauer polynomials in £ = tanh%: 



cm 



(-l)^TT 



iT(a+hn) 



n)\T(a) 



iF 1 



— ~n, ~n + a 



lf n r(a+|(n+l 



n even. 



ifn-l))!r(a 



(61) 



-|(n-l), ^(n+l)+a 
3 
2 



n odd. 



Then, for a 



n + \ 



oo, the hypergeometric polynomials simplify: for 
1 "• 1 • " ' 2x — H n (y/a£), and similarly for 
n odd. Replacing this into Eq. (|59|), with £ = tanhx « sinhx — X\/R and 



s 

n even, 2 -Fi(-|n, a; A; £ 2 ) « i-Fi(-±n; ±; a£ 
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cosh s+n x ~ exp(— s tanh 2 \x)i we obtain the harmonic oscillator wavefunc- 
tions on flat space, 

l j=r n {x) ~ I 1 e-^?/ 2 H^y/iEJxi). (62) 

^ W2™n! ^vr/ziu; 

The factor V^R restores the proper normalization on the x\ axis. 

In addition to the bound states, the sech-trough Poschl- Teller potential 
also has free states with energy above the asymptotic value of the poten- 
tial lim^-too V(x) = \\iio 1 R 2 . These scattering solutions contain associated 
Legendre polynomials of imaginary upper index: 

= 2n P ^ tanh ^ a = \± Ju?/# + 1. (63) 
These wavefunctions are Dirac-orthonormal. 



4.2 Momentum representation of the wavefunctions 

The bound wavefunctions of the Poschl- Teller sech-trough potential in the 
momentum representation, ip^ip) are f° un d through the ordinary Fourier 
transform [Eqs. (|11f)-(|T9|) for D = 1 and ([)!])] of the wavefunctions V'n(x) 
found in Eq. (|5S|). The result is 



R f°° 

— / dxexp(-ipRx)ipn(x) 

Z7T J-oo 



T(2s -n + 1) l r (|( s ~ n ~ l P R ) 



R 

2 \ 7r(s — n)n\ T(s — n) 2 

—n, 2s — n+1, |(s — n — ipR) 
s — n + 1, s — n 



X3F2 



i) n R yJ(s-n)n\T(2s-n+l) 



(64) 



r(|(s — n — ipR 



2^ r(s)r(s + i) 

xR n (-±ipR; |(s— n), n), |(s— n), |(s— n)+l), (65) 
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where Rn(z; a, /3, 7, 8) are the continuous Hahn polynomials [pTO | . On the 



other hand, the unbound solutions (p3| ) are not square-integrable, so their 
Fourier transform must be performed allowing for the phase difference be- 
tween asymptotic incoming and outgoing waves, as determined by the scat- 
tering properties of this Poschl- Teller potential. We shall not further detail 
the free states here; they can be found in Ref. || among the coupling co- 
efficients between the D + x D~ — > + / C irreducible representations 
series of SO(2,l). 

4.3 Wigner function for the oscillator eigenstates on 
the hyperbola 

On the one-dimensional hyperbola 7i\, the Wigner function (ES) collapses 
to (H)-®, its usual form in Eqs. ©-© for D = 1. 

We are interested in the single-function form W(f\x, p) = W(f, /|x, p) = 
W(f , f\p, —x) for the Poschl- Teller wavef unctions, whose explicit form is 
in Eq. ([^) for i[>n(x), and in Eq. (f>|) for 4>n(p); we find the latter more 
amenable to analytic solution. We change the integration to a contour along 
the imaginary axis, and find 

rxrf is\ \ r2 s ~ n x^ (-n) m (-n)i 



8tt 2 n\ r(2s-n+l) n fji Q T(s-n+m) T(s-n+l) 

r(2s-w+m+l)r(2s-7i+Z+l) 1 
X T(s-n+m+l)T(s-n+l+l) l\m\ MPh 1 j 

where i^(x>p) = — — / dzY{\{s—n)—\ipR—z) T(|(s— n)+^ipR+z+m) 

11 J —too 



X 



e~ 4xz T(hs-n)+kpR-z) Y(\(s-n)-\ipR+z+l). (67) 



The last integral can be computed in the complex plane straightforwardly, 
and leads to a pair of complex conjugate 2 -F\-functions of e~ 4x . We thus 
finally have 
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w{r n \x,p) 



2R (s-n) n\ e' 2 ^-^ 
~ irT(2s-n+l) 

" T(2s-n + m + l)T(2s-n + l + l) 
X m ^o r(s-n+Z + 1) r(s-n+m+l) T(s - n + Z) 

Z! m! (n — m)! (n — /)! ^ ^ 

x Re J r(ipi2) T{s-n+l-ipR) e 2ipRx 



x 2^1 



' s — n + m, s — n + I — ipR, 
1 — ipR 




We note that this expression is suitable for numerical computation only for 
X > because it converges fast, but it holds everywhere analytically, with 
the reflection symmetries W \^\x, p) = W(ip^\— XiP) — W(ip^\x, —p)- 

The Wigner functions, together with their marginal projections I'i/r^)! 2 
and \ipn(p)\ 2 i are shown in Figures 1, for n — 0, 1, 2, 3, and for the potential 
depth parameter s = 4 and 30. It can be appreciated that the Wigner 
function of the most tightly bound states resemble the familiar Gaussian- 
bell form of the harmonic oscillator ground state. According to (|60"1), for 
s = 4 there are only 5 bound states (n = 0, . . . , 4), and as the energy of the 
state approaches the binding energy, the wavefunction stretches in space with 
ever-smaller momentum in a neighborhood of the classical turning point. The 
contraction limit can be appreciated in the s = 30 column, corresponding 
to a large binding energy; the Wigner functions for the eigenstates in Eq. 



(59) approach the familiar Laguerre- Gaussian form that corresponds to the 



Wigner function of the harmonic oscillator wavefunctions on flat space. 



5 Concluding remarks 

We have generalized the Wigner quasiprobability distribution function by 
replacing the oscillating exponential functions of the standard version, which 
are Dirac solutions of the free Schrodinger equation on flat space, by the 
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Shapiro functions, because they are solutions to the Laplace-Beltrami equa- 
tion on a simply connected hyperbolic space, while respecting the midpoint 
condition through an appropriate Dirac-like A D (x;x',x"). The role of the 
Fourier transform in the standard version is transfered by the Shapiro func- 
tions to a Mellin-like transform between the position and momentum coordi- 
nates of phase space. Indeed, the relation between what were called position 
and momentum variables is actually defined by the argument and index of 
the basis of Shapiro functions. Thus built, the proposed Wigner function is 
covariant under the group SO(D, 1) of motions of the hyperbola, with the 
hyperbolic translations extracting a multiplier factor. The correct marginals 
are found and the contraction limit to flat space returns the standard Wigner 
function. 

The transformations of the momentum direction under translations of the 
hyperbola are the (unique) action of SO(D, 1) on the sphere S D [[TTJ. There 
are several models of Hamiltonian systems where momentum is restricted to 
a sphere, such as geometric and Helmholtz (monochromatic) optics [ 171, HI]- 
In the first model, a 'Descartes sphere' of momentum vectors p which is of 
radius |p| = n(x) (the refractive index), is associated to each point x in space. 
This sphere of momentum vectors has been subject to Bargmann's 'boost' 
transformation in Refs. |19| (which is unique for Lorentz groups acting on 
spheres) and here given by fl43|), with a canonically conjugate transforma- 
tion of the ray positions at a plane screen. The resulting phenomenon has 
been called relativistic coma, because the aberration in the images projected 
on moving screens in vacuum is comatic. The Doppler effect is of course 
absent, so this transformation cannot be called physically relativistic in the 
Helmholtz case, which involves a configuration space with a nonlocal metric. 
Indeed, the difference between the various models consists in their space of 
positions; here, it is the hyperboloid. A Wigner function on spheres will be 
examined in part II of this title; it is expected to clarify further the use of 
function bases to define conjugate variables as a simile or substitute for phase 
space. 

We point out that the example of the one- dimensional Wigner function of 
Poschl- Teller wavefunctions is of interest not only for the traditional quantum 
mechanical model, but also for the paraxial propagation of light along shal- 
low non-harmonic waveguides whose index of refraction has a sech 2 profile, 
as given by Eq. fl22"|) . Finally, one of the manifestations of higher symme- 
try is the appearance of 'closed-form' wavefunctions expressed in terms of 



20 



well-known (and some not-so-well known) special functions. Generally, the 
Fourier transforms and Wigner functions of such wavefunctions are again 
known special functions because symmetry, when it occurs, is displayed in 
phase space. 
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Figure caption 



Figure 1. Wigner functions of the Poschl- Teller eigenfunctions Vn(x) 011 a 
quadrant of phase space (axes are position X\^s and momentum pR/ y/s; the 
quadrants have reflection symmetry across the axes). Rows are numbered 
by the mode n — 0, 1, 2, 3. Left: s = 4 (so only five states, n — 0, . . . , 4, are 
bound); right: s = 30 (so states are bound up to n = 30). White is the maxi- 
mum, black is the minimum; the shade at the upper right corner corresponds 
to zero. From each Wigner function we project up the marginal distribution 
of position |^(x)| 2 , and right the marginal of momentum \ijj^(p)\ 2 . 
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